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ABSTRACT 



Baryon Acoustic Oscillations (BAO) provide a "standard ruler" of known 
'"5 ■ physical length, making it one of the most promising probes of the nature of dark 

' ' energy. The detection of BAO as an excess of power in the galaxy distribution at a 

?-j ■ certain scale requires measuring galaxy positions and redshifts. "Transversal" (or 

"angular") BAO measure the angular size of this scale projected in the sky and 
provide information about the angular distance. "Line-of-sight" (or "radial") 
'nI" ■ BAO require very precise redshifts, but provide a direct measurement of the 

ly^ ■ Hubble parameter at different redshifts, a more sensitive probe of dark energy. 

f^ ' The main goal of this paper is to show that it is possible to obtain photometric 

Q ■ redshifts with enough precision (cr^) to measure BAO along the line of sight. 

t~^ ■ There is a fundamental limitation as to how much one can improve the BAO 
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measurement by reducing cTz- We show that Oz ~ 0.003(1 + z) is sufficient: a 
much better precision will produce an oversampling of the BAO peak without 
a significant improvement on its detection, while a much worse precision will 
result in the effective loss of the radial information. This precision in redshift 
can be achieved for bright, red galaxies, featuring a prominent 4000 A break, 
by using a filter system comprising about 40 filters, each with a width close to 
100 A, covering the wavelength range from ~ 4000 A to ~ 8000 A, supplemented 
by two broad-band filters similar to the SDSS u and z bands. We describe 
the practical implementation of this idea, a new galaxy survey project, PAltj, 
to be carried out with a telescope/camera combination with an etendue about 
20 m^ deg^, equivalent to a 2 m telescope equipped with a 6 deg^-FoV camera, 
and covering 8000 sq. deg. in the sky in four years. We expect to measure 
positions and redshifts for over 14 million red, early-type galaxies with L > L^ 
and iAB ^ 22.5 in the redshift interval 0.1 < 2; < 0.9, with a precision az < 
0.003(1 + z). This population has a number density n > 10"'^ Mpc~^h^ galaxies 
within the 9 (Gpc/h) volume to be sampled by our survey, ensuring that the 
error in the determination of the BAO scale is not limited by shot-noise. By itself, 
such a survey will deliver precisions of order 5% in the dark-energy equation of 
state parameter w, if assumed constant, and can determine its time derivative 
when combined with future CMB measurements. In addition, PAU will yield 
high-quality redshift and low-resolution spectroscopy for hundreds of millions of 
other galaxies, including a very significant high-redshift population. The data set 
produced by this survey will have a unique legacy value, allowing a wide range 
of astrophysical studies. 

Subject headings: large-scale structure of universe — cosmological parameters 



1. Introduction 

Physical Cosmology has recently entered the precision era. This transition has been 
propelled by the gathering, over the past decade, of unprecedented high-precision data sets 
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for several cosmological observab les. The combined an a lysis of the Cosmic Mi c rowave Back- 



groun d (CMB) anisotropies (e . g..lSpergel et al 



fl2006[ ): lRedhead et al. lfl2004} ): lReichardt et al. I (120081 ): [Dickinson et al 



2007f ): lKomatsu et al. I (120081) : I Jones et al. 



scale measurements at incr easingly higher redshifts (e.g.. 



2004 



Percival et al. 



Wood-Vasev et a 



(120061) : iRiess et al. I( 20071) ) andprobesof large-scale structure (jCole et al. 



(I2004f)') -w i th distance- 
(l2007h:lAstieL^laL 



2005l : lTegmark et al. 



20071 : IXegmark et al. I l2006l : iHiitsi I l2006al ) yields a remarkably consis 



tent picture: a spatially flat universe that has started a phase of acceleration of the expansion 
at the present epoch. From the observations gathered so far this acceleration is consistent 
with the effect of a cosmological constant, but it may also be caused by the presence of 
a dynamical energy component with negative pressure, now termed Dark Energy (DE), or 
might also point to a fundamental modification of our description of gravity. The answer to 
what is the exact cause is likely to have profound implications for cosmology and for particle 
physics. 



Two recent collective reports, one by the US Dark Energy Task Force (DETF) (jAlbrecht et al 



20061 ) . convened by NASA, NSF and DOE, and another by the European ESA-ESO Working 
Group on Fundamental Cosmology (jPeacock et al. II2006I ). have identified the most promis- 
ing observational strategies to characterize DE properties in the near future. These reports 
concluded that the method based on measurements of Baryon Acoustic Oscillations (BAO) 
from galaxy redshift surveys is less likely to be limited by systematic uncertainties than other 
methods that are proposed. It appears that, while recognizing the need for a combined strat- 
egy involving two or more independent techniques, BAO measurements can substantially 
contribute to increase the accuracy on the DE equation of state. 

Baryon Acoustic Oscillations are produce d by acoustic waves in the photon-baryon 
plasma generated by primordial perturbations (JEisenstein fc Hu Ill998l ). At recombination 
{z ~ 1100), the photons decouple from the baryons and start to free stream, whereas the 
pressure waves stall. As a result, baryons accumulate at a fixed distance from the original 
overdensity. This distance is equal to the sound horizon length at the decoupling time, 
tbao- The result is a peak in the mass correlation function at the corresponding scale. First 
detections of this excess were rec ently reported, at a significance of about th r ee standard- 



deviations, both i n spectroscopic (JEisenstein et al. 



and photometric (jPadmanabhan et al. 1 120071 : iBlake et al. 1 120071 ) galaxy redshift surveys 



2005 



Percival et al. Il2007bl : iHiitsi II2006I ) 



The comoving BAO scale is accurate ly determined by CM B observations {tbao = 
146.8 ±1.8 Mpc for a flat AC DM Universe (IHinshaw et al. Il2008[ l). and constitutes a "stan 
dard ruler" of known physical length. The existence of this natural standard ruler, measur 
able at different redshifts, makes it possible to probe the expansion histo ry of the Universe, 



and thereby dark energy properties and the Universe geometry (see, e.g.. lSeo &: Eisenstein 
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fl2003h : iBlake fc Glazebroo 



forts t o measure BAO ( e. 



^ 



{ 20031) and re f erences there i n). This motivate s the present ef- 



ADEPT 



f 2008h.lDESlfl2008 



fl2008h . IsPACeI J2008h . Isasset et aL I J2005h . IWigdeZ I J2008h l 



HETDEX I j2008h . |Pan-STARRS 



Broad-band photometric galaxy surveys can measure the angular scale of tbao in several 
redshift shells, thereby determining dA{z)/rBAo, where dA^z) is the angular distance to the 
shell at redshift z. If galaxy redshifts can be determined precisely enough, the BAO scale can 
also be measured along the line of sight, providing a direct measurement of the instantaneous 
expansion rate, the Hubble parameter (or actually of H{z) vbao), at different redshifts. This 
quantity is more sensitive to the matter-energy contents of the universe compared to the 
integrated quantity dA^z). The direct determination oi H{z) distinguishes the BAO method 
from other methods. In addition, since systematic errors affect the radial and tangential 
measurements in different ways, the consistency between the measured values of H{z) and 
dA^z) offers a test of the results. 

As a rule of thumb, in order to get the same sensitivity to the dark-energy parameters, 
a galaxy redshift survey capable of exploiting the information along the line of sight needs 
to cover only ~ 10% of the volur ne covered by a comp arable survey that detects the scale 
in the transverse direction only (JBlake &: Bridle II2005I ). When covering a similar volume, 
precise enough redshift measurements can provide substantially tighter constraints on the 
DE parameters. 

Large volumes have to be surveyed in order to reach the statistical accuracy needed to 
obtain relevant constraints on dark-energy parameters. Enough galaxies must be observed 
to reduce the shot noise well below the irreducible component due to sampling variance (see 
section [27il) . The usefulness of the correlation along the line of sight favors spectroscopic 
redshift surveys that obtain very accurate redshifts, but the need for a large volume favors 
photometric redshifts that can reach down to fainter galaxies. 

The intrinsic comoving width of t he peak in the mass correlation function is about 
15 Mpc/h, due mostly to Silk damping (jSilk Ill968l ). This sets a requirement for the redshift 
error of order a{z) = 0.003(1 + z), corresponding to 15 Mpc/h along the line of sight 
at 2; = 0.5. A much better precision will result in oversampling of the peak without a 
substantial improvement on its detection, while wor se precision will, of cou rse, result in the 
effective loss of the information in the radial modes (jSeo fc Eisenstein 1120031 1. Note also that 
in the presence of substantial redshift errors, the error distribution needs to be known and 
accurately corrected for when inferring the BAO scale. 

It has usually been assumed in the lite rature that photometric redshifts are not pre - 
cise enough to measure "line-of-sight" BAO (ISeo fc Eisenstein II2003I : iBlake &: Bridle II2005I ). 
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While this is true for broad-band photometric surveys, here we examine how one can reach 
the required redshift precision with narrow-band photometry. We find (see section [3]) that 
redshifts of luminous red galaxies can be measured with a precision a{z) ~ 0.003(1 + 2) using 
a photometric system of 40 filters of ~ 100 A, continuously covering the spectral range from 
~ 4000 to ~ 8000 A, plus two additional broad-band filters similar to the u and z bands. 

We describe the practical implementation of this idea, a photometric galaxy redshift 
survey called PAU (Physics of the Accelerating Universe). PAU will measure positions and 
redshifts for over 14 million luminous red galaxies over 8000 deg^ in the sky, in the range 
0.1 < z < 0.9 (comprising a volume of 9 (Gpc/h)^), and with an expected photometric 
redshift precision a{z) < 0.003(1 + z). This redshift precision makes it possible to measure 
radial BAO with minimal loss of information. The PAU survey can be carried out in a four 
year observing program at a dedicated telescope with an effective etendue ~ 20 m^deg^. 

The outline of the paper is as follows. Section [2] discusses the scientific requirements for 
such a survey. An optimization of the survey parameters follows in section [3l while section H] 
presents a possible baseline design for the instrument. The science capabilities of the survey 
are given in section O Finally, section E] contains the summary and conclusions. 



2. Scientific Requirements 

Given current priors on other cosmological parameters a measurement of the expansion 
rate history H{z) with percent precision will translate into a measurement of the DE equation 
of state w of few times this precision. We will focus here on how well we can do this by 
measuring the BAO feature at comoving size tbao — lOOMpc/h and use it as a standard 
ruler in both the radial and tangential directions. 

We first explore how an error in the BAO scale vbao translates into an error in the 
DE equation of state w. The relation is different if we measure the scale in the radial or 
in the tangential (ie angular) directions. We then move to show how the BAO scale can be 
measured statistically using galaxy surveys. We start with a visual ilustration of the problem 
and a brief presentation of two N-body simulations that we used in order to study the main 
goal of this paper. We then show how this scale can be measured using the statistics of galaxy 
density fluctuations and relate the error in tbao to the volume of the survey, given perfect 
distance indicators. We flnish by considering what redshift precision is required to maintain 
a given precision in the vbao measurement and how this is limited by different systematic 
effects. Note that we focus here in showing the implication of the photo-z accuracy in the 
measurements of the monopole, which combines the radial and perpendicular information. 
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These two components can be separated by considering the anisotropic correlation function 
(see Okumara et al. 2007, Padmanabhan & White 2008, Cabre & Gaztanaga 2008, and 
Figs. 17-18 in Gaztanaga, Cabre & Hui 2008, which show how the anisotropic correlation 
function changes for PAU-like photo-z precisions). 

We note that these are just rough estimates to show the viability of this approach. In a 
real survey, there might be other sources of systematic errors that have not been taken into 
account in detail here. However, experience indicate that the actual data themselves can be 
used to study and minimize those sources of errors. 

In what follows we as sume a flat FRW cosm ology, with cosmological parameters com- 



patible with WMAP data flHinshaw et al. 1120081 ). 



2.1. BAO scale and DE equation of state 

In a galaxy survey we measure distances in terms of angles and redshifts. These observed 
quantities can be related to known distances, such as rsAo, using the FRW metric. The 
differential radial (comoving) distance is inversely proportional to the expansion rate H{z) = 
d/a: 

dr{z) = J— dz , (1) 

while the angular diameter distance is proportional to the integral of dr{z): 

for a fiat universe. In particular, measurements of the charateristic size of the BAO feature 
in the radial (Szbao) and tangential {S9bao) direction relate to the known comoving BAO 
scale vbao as: 

ozbao = rBAO [o) 

c 

dA{z){l + z) 

We therefore have, neglecting for a moment the uncertainty on the determination of vbao by 
CMB observations, that a relative error in the measured size of the BAO feature corresponds 
to a relative error in either dA{z) or H{z): 

J, _ a{60BAo) _ ^dA 



^BAO 



59bao dA 



.L _ (rjSzBAo) _ Aij" 

^BAO = —FZ - -JT [^) 

OZbaO -n 
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For this argument, we will write the expansion rate as: 

H\z)/Hl = ^m{l + zf + (1 - fi„) (1 + zf^^+^^ (6) 

which corresponds to a fiat universe (f2m + ^de = 1) with a constant equation of state 
w = p/p, with p the pressure, and p the density of the dark energy. Figure [1] shows the 
relative change in Eq.Q (in percent) as a function of the relative changes in w with respect 
to w = —1 for different redshifts. We show the cases for z = 0.3 — 1.0 which will be relevant 
for our study. As can be seen in the figure, a 1% error in Abao (our goal) at 2; ~ 1.0 results 
in a ~ 4% uncertainty in w, but the precision varies with redshift. A constant 1% error in 
the angular distance quickly degrades the w precision with decreasing redshift, from 4.5% 
aX z = 1.0 to 8.0% at z = 0.3, while the radial distance achieves a more uniform precision in 
w, 3.5 — 4.5% in the whole redshift range. 

This illustrates the advantage of having a good radial measurement. Angular distances 
provide a good geometrical test, while radial distances tell us directly the instantaneous 
expansion rate. In addition, comparing relative sizes of the BAO feature when measured 
parallel and perpendicular to the hne of sight will provide us with a consistency test. 



2.2. BAO scale in N-body simulations 

To support some of the main claims of this paper we have used large N-body dark matter 
(DM) simulations, using the MICE collaboratioiu set-up. In particular, we have computed 
non-linear DM clustering statistics in terms of the 2-point correlation function, ^(r), and its 
Fourier transform, the power spectrum, P{k), and we have assessed the impact of real world 
systematic effects on the BAO mesurements based on these standard estimators. 



MICE simulations have been run using the Gadget-2 code (jSpringel 1 12005| ) on the 



Mare Nostrum supercompu ter at BSC[J||, with a modification to produce outputs in the light- 



cone (iFosalba et al. 1120071 ). We focus here on two simulations, shown in Table [H correspond- 
ing to a fiat concordance ACDM model with fi^ = 0.25, ^a = 0.75, fi;, = 0.044, n^ = 0.95, 
(Ts = 0.8 and h = 0.7. Halos were obtained from the z = 0.5 comoving output using the 
Friends-of- Friends (FoF) algorithm with linking length 0.164. The larger simulation has a 
dynamic range close to five orders of magnitude. 

Figure [2] shows a thin slice of the light-cone built from the MICE3072 simulation. We 



^http://www. ice. cat/mice 



"Barcelona Supercomputer Center, http://www.bsc.es 



build the light-cone placing the observer at the origin, so that cosmic time (redshift) cor- 
responds to the radial direction, which expands from z = to z = 1 (corresponding to a 
comoving radius of 2400 Mpc/h for our cosmology). The bottom panel corresponds to the 
true dark matter distribution in real space. The next panel up shows the redshift space 
distribution, where the radial positions are distorted due to peculiar (gravitationally in- 
duced) motions away from the Hubble flow. In order to model this distortion, we add the 
radial component of the peculiar velocity Vr of each particle to its (real space) position: 
s = r + fvr{l + z)/H{z), where / ~ 1 for the assumed ACDM cosmology at 2 ~ 0.5. We 
note that, in this image, distortions can only be seen whenever they are much larger than 
the pixel size Ar^j^. ~ 3Mpc/h or, in velocity units, Avpix — 300km/s. This implies that the 
so called Fingers of God effect (see explanation below) cannot be detecte d because it a rises 



from random peculiar velocities of order Avpix- Instead, the Kaiser effect (iKaiser 1119871 ) due 
to a large-scale coherent infall is visible as an enhancement of the filamentary structures 
perpendicular to the line of sight. 

The two top panels include in addition a radial distortion due to photo-z errors which 
we assumed to be Gaussian distributed. Thus, they can be modeled by randomly displacing 
the particles along the line of sight according to the probability distribution, 

f{6r^) ~ exp [-(1/2) {6rJA^f] , (7) 

with the smoothing scale A^ and the photo-z error a^ related through the Hubble parameter 
H{z) as Az = az{l + z)c/H{z). Third panel up assumes az = 0.003(1 + z), which is roughly 
the photo-z error expected for PAU galaxies, while the top panel corresponds to an order- 
of-magnitude worse case, a^ = 0.03(1 + z). 

Overall this figure illustrates how the image is degraded both by peculiar velocities and 
redshift errors. Comparing the middle panels, it is evident that redshift errors produce, on 
average, much stronger distortions than peculiar velocities. 



2.3. Redshift errors 

We now turn to a more quantitative estimate of the minimum radial resolution required 
to detect the BAO scale in 3D. 

The BAO signature appears in the two-point correlation of particles as a single bump 
at a scale tbao — 100 h~^ Mpc with an intrinsic width ArsAo — 10 h~^ Mpc and relative 
amplitude of about a factor two with respect to a non-BAO model with the same broad-band 
shape. Therefore, a simple approach would be to look at the degradation of this peak as we 
increase the redshift error. 
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To be more realistic, we have studied the clustering of halos (selected with a FoF al- 
gorithm with linking length of 0.164) since they are more closely related to the observed 
galaxies and clusters of galaxies than DM particles. In particular, Luminous Red Galaxies 
(LRGs) are thought to populate large DM halos and closely trace the halo distribution. 

At large scales the halo and matter density fluctuations are related by a linear bias 
factor {6 halo = b6matter), which translates into a 6^ scaling for the 2-point funcion. Halos 
with large mass cut in Tabled] have 6 ~ 3, which was chosen to magnify possible non- linear 
effects. To estimate errors (eg. in Fig. Hj), we use instead a lower mass cut which corresponds 
to 6 ~ 2 ar id better matches the clustering of LRGs tha t have already been used to measure 



BAG fe.g. lEisenstein et al. II2005I : iTegmark et al. 1120061 ) 



Figure [3] shows the two-point correlation function (black empty circles) traced by all 
halos in MICE3072 at 2; = 0.5 (see Tabl^. The dashed line corresponds to the linear 
correlation and the black solid line to the nonlinear model of Renormalized Perturbation 



Theory (RPT) (ICrocce &: Scoccimarro II2008I ). both biased with 6 = 3. As illustrated by 



Fig. El (solid black line), the simple scale independent linear bias scheme works reasonably 
well, but this assumption certainly needs to be tested more accurately. 

We then estimated the impact of photo-z errors using Gaussian distortions in the radial 
direction as described in the previous section. The corresponding correlation function of 
this smeared distribution of halos is shown in Fig. [3] for different values of the photo-z error, 
cr^/(l + z) = 0.003 (as expected for PAU) 0.006,0.03, with red, blue and green symbols 
respectively. 

In turn, this smearing can be modeled in Fourier Space by damping the power spectrum 
along the line-of-sight as 

P,(fc,/i) = b^Pni{k)exp [-k^Ay] , (8) 

where Pni is the nonlinear power spectrum from RPT, the linear bias is b = 3, and fi is 
the cosine of the angle with the line of sight. The red, blue and green solid lines in Fig. [3] 
correspond to the angle averaged Fourier transform of Eq. ([H]) for a^ = 0.003, 0.006 and 0.03 
respectively. 

In summary. Fig. [3] illustrates that one can basically recover the right BAG shape once 
the error is better than about 0.003(1 + z). Larger errors erase the BAG bump and will 
result in the loss of cosmological information. The change can be roughly quantified by the 
ratio between the amplitude of the BAG peak (at r ~ 108 Mpc/h) and the amplitude in 
the valley (at r ~ 85 Mpc/h). For r'^C,{r) in the right panel of Figure [3l this ratio is about 
1.8 in real space with no photometric errors and decreases smoothly to 1.5 as we increase 
the error towards 0.003(1 + z). For larger errors, this ratio decreases more rapidly an gets 
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all the way to unity for 0.006(1 + z), as shown in the figure. This makes sense because 
0.003(1 + z) corresponds to a comoving scale of about 15 h~^ Mpc at 2; = 0.5, which matches 
the intrinsic width (Silk damping) of the BAO peak. Although this is just a rough estimate, 
it is all we need as a starting point for our considerations below, and it agrees with other 
co nsiderations based on count in g the number of modes i n 2D and 3D surveys and the work 



of lBlake fc Glazebrook I (120031 ): ISeo fc Eisenstein I (l2003f ). 



The clustering analysis presented above is in real rather than redshift space. Redshift 
space distortions can be modeled as a combination of two separate effects: coherent an d 



random peculiar velocities. The first term is the so-called Kaiser effect ( iKaiser 1 119871 ). 
which increases the amplitude of clustering at large scales by a factor ~ (1 + /9/i^)^ (where 
(3 = Q'^/b, b is the bias and fi the cosine of the angle with the line of sight). Our analysis 
allows for such an effect by incorporating a larger effective bias in the correlation function 
monopole (e.g. 6 = 3 as in Fig. [3]). The effect of random velocities can be modeled as 
a Gaussian damping, very similar to photo-z errors in Eq.([8]) but where a^ is replaced by 
(Tp/2, where dp is the one dimensional galaxy pairwise velocity dispersion (the factor of 1/2 
is because a velocity difference has twice the variance of a single velocity). The net effect is 
that the density field is convolved with a one dimensional random field with a net dispersion 
a"^ = (7^/2 + cr^. The typical value of Up in our DM simulation is smaller than the photo-z 
errors considered here. This is also the case in regions of high density, populated by LRG 
galaxies, and whe re cxp could be larger, e.g, ap ~ 400km/s/c ~ 0.0013 for r > 5Mpc/h 



( iRoss et al. 1120071 ). This is clearly seen by comparing the two middle panels of Fig 12] where 
it is evident that even in regions of high density the photo-z distortions are larger than 
the redshift space distortions. Thus effectively az > <Jp and we can consider this effect 
subdominant in our considerations. 

We note here that there is a fundamental limitation as to how much one can improve 
the BAO measurement by reducing the photo-z error. The photo-z error a^ — 0.003(1 + z) 
proposed by the PAU survey is close to optimal. Redshift space distortions and non-linear 
effects can produce distortions that are comparable to this value, depending on what is the 
(biased) tracer that is used to measure BAO. 



2.4. Estimating the BAO scale 

Armed with the conclusions from the previous discussion about photo-z errors and other 
nonlinear effects (clustering, bias and redshift distortions), we are now in a position to give 
an estimate of the expected 1-a error determination from a survey with the characteristics 
of PAU. To this end we will employ the 2-point statistics in Fourier Space (i.e. the power 
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spectrum P{k)). 

To estimate the error in the measurement of the power spectrum we will resort to the 
commonly used expression: 



ap 



AP{k) 



Nm{k) 



1 + 



1 



nP{k) 



(9) 



which can be derived from iFeldman. Kaiser fc Peacock I (119941 ) (see, for instance, Martinez 
& Saar 2002, §8.2), where N^ik) is the number of Fourier modes present in a spherical 
shell extending from k io k + Ak. In terms of the survey volume V, we have Nm{k) = 
V {4tt k"^ Ak) / {27r)^ . The first term in Eq. iQ corresponds to the sampling error and is inde- 
pendent of redshift. The second term corresponds to Poisson shot-noise, and n denotes the 
number density of observed galaxies in the survey. This formula is exact when the probabil- 
ity density function of spectral ampli tudes is Gaussian an d a very good approximation when 
the shot-noise term is negligible (see lAngulo et al. II2008I ). 



For k < 0.12/i^^Mpc and for halos that host LRGs we expect that P{k) > 2 x 
10^/i~^Mpc (6 ~ 2) at z ~ 0. This agrees well with actual mea surements of P{k) for 
LRGs in the SDSS catalogue (see e.g. Fig. 4 in lTegmark et al. I (120061 )). As we will show be- 
low, the PAU number density of LRGs is expected to be n > 0.001 h^ Mpc~'^ for z < 0.9 (see 
FigHH below), which implies that the Poisson shot-noise contribution to the error in Eq.Q 
is smaller than 8% at z ~ 0.5, even taking into account the degradation from a photo-z error 
of 0"^ = 0.003 as discussed before. 

The bump in the spatial correlation function translates into the power spectrum as a 
series of damped oscillations of a few percent in relative amplitude. This is shown in Fig. |4] 
that contains the power spectrum of DM (left panel) and halos of mass M > 4.7 x 10^^ h^^ Mq 
(right panel) measured in the MICE1536 simulation. In both panels the measured spectra 
have been divide d by a smoothed one with the same broad band power obtained from the 



data themselves (jPercival et al. Il2007bl : lAngulo et al. II2008I ). The solid red lines in Fig. H] 



show that this ratio can be roughly modeled as, 

P{k) ^ 1 + Akexp[-{k/0.1hMpc-^)'^]sm{rBAok), 
which illustrates how P{k) depends on the BAG scale. 



(10) 



Moreover, the discussion in Sec. 12.31 that led to Fig. [3] validates to a good extent that 
photo-z, clustering, bias and redshift distortions can be modelled in the power spectrum 
monopole as the angle average of. 



P{k,f,) = b'PUm + Pfi'Yexp [-fc^AV] 



(11) 
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where A = cr(l + z)/H{z) and a = \/cr'^ + crp/2. Therefore, except from nonlinear cluster- 
ing, which is stronger than the intrinsic Silk damping of BAO, and is responsible for the 
exponential damping in Eq. (TTOj) . the remaining effects are multiplicative contributions to 
the measured P{k), and they factor out when constructing P. This means that the above 
described systematic effects do not affect the BAO signal in the spherically averaged P{k), 
although they do increase the associated errors as we will discuss below. Thus Eq. (ITOl) 
allows to compute how the measured P varies with the BAO scale, 

dP/drBAo - Ak"^ exp[-(A;/0.1/iMpc"^)2] cos {rBAok) (12) 

Such variation in P produces a shift in the x^ fitting to measurements of P{ki) given 



by, 



--E^-U(;£,)^(a'^""- 



where, 



, 1 r'- ^<^exp[-2{x/10Mr]cos\x) 

^ ^ " (2vr)2 io (1 + l/nPy ' ^"^ 



m is the number of BAO oscillations included in the fit [k^ax = '^t^^t^/^bao) and ^bao = 
^f BAO /'^ BAO- The shot-noise term includes the full power spectrum given by Eq. ( ITTl) and 
accounts, in particular, for photo-z errors. In deriving Eq. flT3l) we have explicitly used that 
for our reference cosmology tbao = 108.6 h^^Mpc. 

A 1-sigma determination of Abao alone corresponds to Ax^ = 1, so that, 

Abao\ax^=i = -^rf- J, ./., r- (14) 

\ V J I[m\{AlrBAo) 

From Fig. m we find that A/tbao ~ 0.02 fits well both halos and dark matter clustering. We 
then assume m ~ 2.5 which corresponds to the range < /c < 0.14 /i Mpc~^ (including two 
BAO peaks) and obtain from Eq. ([1] 



iBAo — 0.33 



/ 3 \ 1/2 

r-^\ ~0.33%-yi3/i-3GpcVy (15) 



when we neglect completely the shot-noise term in Eq. flTT|) . Including the shot-noise but no 
photo-z degradation yields a prefactor of 0.35% (for ap = 400km/s). If we also add a photo-z 
error of a^ = 0.003 (as PAU) we find 0.36%. If instead, we add a photo-z error of a^ = 0.03, 
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the amplitude rises to 2.2%. In other words, we expect PAU to yield a measurement of 
BAO with only a 10% degradation with respect to an ideal survey, whereas a survey with an 
order-of-magnitude larger photo-z is expected to produce a factor ~ 6.5 worse measurement. 
According to this, in a 3D analysis, the relative error in the BAO scale is just approximately 
equal to the inverse of the square root of the number of independent regions of size t^ao 
that are sampled by our survey, and it is quite robust in front of close to optimal (PAU-like) 
photo-z error and nonlinear effects. For V ~ 10h~^ Gy cf we get abou t A ^ An — 0.5%. The 
above esti mate is in good agree r aent with Table 2 in lAngulo et al. I (120081 ) and with the 
analysis in iBlake fc Glazebrook I (120031 ) and ISeo fc Eisenstein I (120031 ) . 



If we limit ourselves to optical surveys of LRGs, we have z < 1. To get to ^ ~ 
10 h^^ Gpc^ we will have to map of the order of 8000 sq. deg. There are roughly two million 
LRGs with luminosity L above the ch aracteristic galaxy luminosity L* in 1000 sq. deg. at 
z < 0.9 with magnitude Iab < 22.5 (iBrown et al. 1120071 ). However, not all galaxies in a 
given volume need to be measured as long as nP > 3, so that shot noise is sub-dominant in 
Eq. (j9]). We will show below that it is in fact possible to get to nP > 10 with the subsample 
of PAU LRGs that have good quality redshifts. 



3. Survey Simulations 

The main distinctive feature of our survey is the use of photometric information to 
achieve the highly accurate redshift measurements needed to characterize the line-of-sight 
BAO signature. Since such an observational program has not been attempted before, we need 
to prove, at least conceptually, that it is possible to achieve precisions of ^"^/(l -|- 2;) ~ 0.003 
with photometric data. According to section [2] and particularly Eq. 1^, we only need to 
reach such a precision for a galaxy population tracer with a space number density which 
satisfies nP{k) > 3. Ideally we would like galaxies that are luminous so we can observe them 
up to high redshift and that present a spectral energy distribution with distinctive features to 
achieve accurate photometric redshifts. The most luminous of the early type galaxies (LRGs 
or Luminous Red Galaxies) constitute such a population. Their space number density is 
high enough. They are highly biased and they feature a prominent 4000 A break in their 
spectrum which, together with other spectral features, makes possible a precise estimation 
of their redshifts using photometric measurements (Fig. E]). As a matter of fact, bright early 
type galaxies we re the subject of the first attempt to estimate photometric redshifts in the 
seminal paper of iBaum I (119621 ). 



Hickson et al. I (119941 ) was the first to propose using intermedi ate band filters as a vi- 



able alternative to traditional spectroscopy. The COMBO-17 survey (IWolf et al. Il200ll . l2003l ) 
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put this idea into practice, using a combination of traditional broad band and medium band 
filter s. COMBO-17 reaches an accuracy of a^ ~ 0.02(1 + z) for the general galaxy popula- 
tion (JHildebrandt et al. II2008I ). and has reached a scatter of a^ ~ 0.0063 for the bright ellip- 
ticals in the Abell 901/902 superclusters. Taking into account the velocity disp ersion of the 



cluster the authors infer an intrinsic photometric accuracy close to 0.004(1 + z) (jWolf et al. 
2003h . 



Benitez et al. I (120081 ) has shown that the most effective way of reaching high photo- 



z precisions is using a system of constant-widt h, contiguous, non-overlapping filters. The 
ALHAMBRA survey (JMoles et al. 1 12005| . |2008| ) has implemented such a filter system, and 
preliminary results for that survey show that it is possible to get close to 0.01(1 + z) photo-z 
accuracy for the general galaxy population. LRGs usually have higher photo-z precisions 
than the rest of the galaxies, (as it happen s with the COMBO-17 data or the SDSS LRGs 
(jOvaizu et al.ll2008l : lD'Abrusco et al. IboOTl )). and it is expected that LRGs in the ALHAM- 
BRA survey will have photo-z errors substantially below cr^ ~ 0.01(1 -|- z). 

In view of these results, it seems reasonable to suggest that a precision a few times 
smaller that 0.01(1 -|- z) can be reached for the LRG population with filters that are three 
times narrower than those of ALHAMBRA, and about 2-3 times narrower than the medium 
band filters in COMBO-17 (which do not fill the optical range contiguously). In what 
follows we try to demonstrate that the observing program required by PAU is feasible, and 
can deliver redshift values with (^^/(l -|- 2;) ~ 0.003 for LRGs. 

In order to qualitatively understand the relationship between measurement uncertainties 
and the accuracy of the redshift estimation we can use a toy, step-like spectrum, fiat in 
wavelength except for a jump by a factor D in the amplitude at 4000 A: i.e. we assume that 
the spectrum has a fiux of F redwards of the break and F/D bluewards o f it. This roughly 
approximates a low resolution version of an LRG ( Eisenstein et al. II2003I ). If we use a set 
of constant width, contiguous filters of width AA, the fiux in the filter that spans the break 
will be equal to 

f = a^ + il-a)F (16) 

where a = {\b — Ao)/AA -|- 1/2 = Rz + k, here A^ = 4000(1 -|- z) is the observed wavelength 
of the break, Aq is the central wavelength of the filter, the "local resolution" is i? = 4000/ AA 
and k = (4000 - Ao)/AA -1-1/2. We then have that 



f = F[iD-'-l)iRz + k) + l] 



(17) 



and 



1 
R 



i-^ 



D 



1-D 
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The error in the redshift roughly depends on the flux measurement error o"/ as 

D Of . , 

^-f ~ im^^-F (19) 



where we have considered that the error in the determination of F is much smaller than (jf. 

Apart from the photometric error, another source of uncertainty in the redshift esti- 
mation is the intrinsic variabili ty of galaxy spectra around its average, e ven within such a 



homogeneous class as LRGs. (jCool et al. I l2006l : lEisenstein et al. II2003I ). We can include 



this in our toy model as an uncertainty in the 4000 A break amplitude D. Using the above 
formulae we get that 



1 
R{lT- 1)2 

The total uncertainty predicted by the toy model is thus, 



(Tzu < „,^ ,,, ^D . (20) 



-^- R^D-lY ^-l+'^^-flFy- (21) 

We can estimate the intrinsic scatter m D to be o"/) = 0.1, as shown below. To check 
the validity of this formula we can use the photometric observations of LRGs with measured 
spectroscopic redshifts. We have downloaded a LRG catalog with spectroscopic redshifts 
from the SDSS website and estimated their photometric redshifts using the LRG template 
described below, measuring an average error of Oz ~ 0.02(1 + z). Eq. fl2Tl) clearly over- 
estimates the error, since at e.g. z = 0.3, where galaxies have typically (Tf/F = 0.05, it 
would predict (using D = 1.8) an error from the template variability of az^^ = 0.044, a 
photometric error of o"^, = 0.03 and a total error of a^ = 0.05, about twice as large as the 
real result a^ = 0.026. This is not surprising, since real galaxies have many features which 
contain redshift information apart from the 4000 A break. Therefore, although Eq. fl2T]) can 
be useful to qualitatively understand the effects of intrinsic scatter and photometric noise 
on the photometric redshift accuracy, it clearly underestimates the precisions which can be 
achieved in practice. 

The application of Eq. [21] to our set up, with <7f/F = 0.1 gives a^ = 0.006, which again 
is twice what we expect. In what follows we will perform a detailed simulation to show 
that it is feasible to reach the photometric redshift accuracy required for our experiment 
(0.003(1 -|- z)) under realistic observing conditions, taking into account the shape of real 
galaxies, the behaviour of the sky background as a function of wavelength and lunar phase, 
and the expected throughput and efficiency of astronomical instruments. 
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3.1. Observational setup and S/N considerations 

In order to simulate the characteristics of the astronomical site where the PAU observa- 
tions will be carried out, we assume that the sky bri ghtness for th e dark phase of the lunar 



cycle is similar to that of Paranal, as measured by iPatat I (2004J ). For the middle of the 
Moon cycle, or "gray" time we use the values of IWalker I ( 19871 ) for Cerro Tololo. Figure [H] 



shows the assumed sky brightness in the standard UBVRI broad bands. However, due to 
the narrowness of our filters, it is necessary to have a good representation o f the small scal e 



structure of the sky spectrum, for this we use the model optical spectrum of iPuxley I (120081 ) . 
the same used for the Gemini exposure time calculator. 

We have written an exposure time calculator for this task. To calculate the full through- 
put of the system we use the La Palma atm osphere at 1.2 airmass, two aluminum reflections 



and the LBNL CCDs (iHoUand et al. II2003I ) quantum efficiency curves. We also approximate 
the throughput of the filter system using the values of the BARR filters produced for the 
ALHAMBRA Survey. The final result is shown in Fig. O We match our results to those 
of the ING exposure time calculator, SIGNAL, using the same observational setup. To 
reproduce their results, which have been calibrated empirically, we have to degrade our the- 
oretical estimates by 25% — 10% (which are basically the values of the empirical corrections 
they use). We have checked our predictions with preliminary results from the ALHAMBRA 
Survey observations and they agree within 10%. We have compared the predictions of our 
simulator with those of DIET o, the Direct Imaging Exposure Time calculator. DIET esti- 
mates 5o" point source limiting magnitudes oi g = 25.74 (dark time) and r = 24.94, i = 24.49 
(grey time) for 442 s exposures, within ^ 1 arcsec^ apertures and 0.8 arcsec seeing. We can 
scale these results taking into account the effective width of our narrow band filters and the 
relative collecting mirror areas (10 m^ for CFHT vs. vr m ^ for our fiducial telescope), cor- 
responding to limiting magnitudes of mF4gs2 = 23.63, mFQ283 = 23.12 and niFi-m = 22.78. 
Our predictions shown in Fig 7 are mF4,982 = 23.77, mF6283 = 22.98 and mF777i = 23.15. 
Most of the differences can be explained by the introduction by DIET of a coefficient which 
attempts to account for the incorrect measurement of the sky background for very small 
apertures and which worsens the S/N by a factor 1.22 for faint objects. 

The simulations below have been carried out assuming that the survey will use a dedi- 
cated 2m-class telescope, with an effective area of vrm^, and a camera with a 6 deg^ FoV. The 
results will remain qualitatively valid as long as the etendue of the final observational setup 
is roughly the same. Most likely the observations will be carried out in drift-scan mode. 
Since there is no need to change instruments, we expect that the observing efficiency will 



http://www.cfht.hawaii.edu/Instruments/Iinaging/Megacain/dietmegacam.html 
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be very high and that only a maximum of two CCD readouts per filt er will be carried out . 



Assuming that the useful time will be similar to that in Calar Alto (ISanchez et al. 1120071 ) 
and that the moonlight will prevent us from taking data during 3 nights per Moon cycle, 
the number of useful hours per year amounts to 1930. Leaving some room for unforeseen 
incidences, we assume that the total number of hours of exposure time per year will amount 
to 1800. For a survey area of 8,000 sq.deg, and with a 6 sq. deg. camera in a period of 4 
years we expect to be able to expose each field a total of 5.4hrs, or 19440s. 

The b est way of measuring accurate colors for photo-z is using relatively small isophotal 



apertures (jBenitez et al. II2004I ) which maximize the S/N of the color measurements, despite 
the fact that such an aperture leaves out a large amount of the flux, and they are therefore 
not optimal for other scientific purposes. In our S/N estimations we assume that we will use 2 
arcsec^ apertures, which enclose about 40% and 64% of the flux of respectively a. z = 0.2 and 
a z = 0.9 L^ galaxy. In the next subsection we explain how we calculate these corrections. 

A crucial question is how to divide the exposure time between the different filters. At 
each redshift, we identify the filter which corresponds to the 4150A rest frame region, and 
try to detect a LRG L^ within a 2 arcsec^ aperture at that redshift with at least a S/N of 
10. We set a minimum exposure time of 120s, and adjust the maximum exposure time in 
each filter so that the total is below 5.4 hours. The resulting exposure times are < 120s for 
filters bluer than F5446, and increase until they reach the maximum exposure time of 861s 
for F7307 and redder filters. The resulting 5a limiting magnitudes are plotted in Fig. [7] and 
Fig. El 



3.2. Intrinsic galaxy variability 

As Eq. f l2T]) shows, it is necessary to understand the intrinsic variability of LRGs in 
order to estimate the photometric redshift accuracy achievable with them. It is well known 
(Eisenstein et al. 2003, Cool et al. 2006 and references therein) that LRG galaxies (with 
L > 2.2L^) are a remarkably homogeneous class. At a fixed redshift, they form a red 
sequence, which varies slowly and regularly with absolute magnitude and environment. LRG 
galaxies in the red sequence present a scatter of only a few percent in the color defined by 
a pair of filters spanning the 4000A break. Therefore, if we know the absolute magnitude 
of a LRG and the richness of its environment, we can predict its broad band colors with a 
precision of at least ag-r ~ 0.035 (Cool et al. 2006). 

It is not clear however, which are the actual variations in the spectral shape of LRGs 
behind this broad band scatter. The SDSS spectrophotometry is not good enough to ac- 
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curately characterize this phenomenon, since its precision (about 0.05 mags in g-r colors, 
according to the SDSS web site and Adelman-McCarthy et al. (2008)) is of the same order 
or even larger than the intrinsic color variation of real galaxies. In addition, the errors in 
the spectrophotometry are bound to be highly correlated and will be much worse at certain 
wavelength regions, like sky lines. The section of the SDSS website which describes the 
quality of the spectrophotometric calibration shows that below 4000A the flux calibration 
error can be as large as 10% . 

Since it is not feasible to use the SDSS spectral information, we have therefore decided 
to use a different approach to characterize this intrinsic spectral variability. Eisenstein et 
al. (2003) split the spectra of LRG into different classes and samples, and looked at the 
differences amongst them using Principal Component Analysis (PCA). They showed that 
most of the variation between these average classes can be ascribed to a single spectral 
component. It is therefore reasonable to assume that the intrinsic variation for galaxies 
of each class, responsible for the red sequence scatter described above, can be modeled 
approximately using the same PCA component. We have therefore generated a mock galaxy 
sample with L > L^ at z = 0.16 using the red sequence described by Cool et al. (2006), and 
the luminosity function described in Brown et al. (2007), and fit their g — r and r — i colors 
using the average template of Eisenstein et al. (2003) and the first PCA component (shown 
in Fig. 8 of Eisenstein's paper). 

The reason to limit ourselves to those two filters is that Eisenstein et al. (2003) only 
provide spectra in the 3650A — 7000A wavelength range, which does not include other filters. 
The comparison with the average Cool et al. 2006 colors show that we have to slightly correct 
Eisenstein's average template to adapt it to the observations, subtracting the first PCA 
component multiplied by 1.74, and that the required variation of the amplitude of the PCA 
component needed to explain the intrinsic scatter around the red sequence is approximately 
1.8 times that necessary to explain the variation of LRGs with redshift, magnitude and 
environment. The average D4000 is 1.81, and the rms around this average value is 0.104. 
Therefore the 4000A break amplitude seems to display an intrinsic scatter of 6% in real 
galaxies. 



3.3. Input early type catalog 

To describe the early type galaxy population we use the luminosity functions described 
in Brown et al. (2007). We populate a 10 sq. deg. area with galaxies following this distribution 
and exclude those which are fainter than L^ (although it is obvious that many of them will 
be detected as well) and fainter than Iab = 23. 
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LRGs are extended objects, and we have to calculate which is the fraction of the total 
flux which falls within our reference 2 sq. arcsec aperture. For this we use the data on 
galaxy sizes and their evolution provided by Brown et al. (2007), assuming, as they do, that 
galaxies can be well represented by a de Vaucouleur profile. The correction ranges from 1.5 
mags at z=0.1 to ~ 0.5 for z > 0.7. The resulting differential numbers counts and redshift 
distribution are plotted with dashed lines in Figs. [9] and [10] (note that we plot the magnitudes 
corresponding to a 2 sq. arcsec aperture). 



3.4. Results 



To simulate our observations we will redshift and integrate under the corresponding 
filter transmissions a spectrum resembling a typical LRG galaxy, and then try to re c over i ts 
redshift using a Bayesian photometric redshift method (BPZ, described in iBenitez I ( l2000l )). 
Obviously in the real world we will not use a single template for all LRGs between < z < 
0.9: their spectra are known to vary with redshift and luminosity (Eisenstein et al. 2003; Cool 
et al. 2006). However, as that paper shows, the variation is smooth and easy to parameterize. 
This is confirmed by HST very deep observations of galaxy cluster s, where the scatter aro und 
the red sequence remains small (~ 0.03) up to z = 1 and higher (JBlakeslee et al. II2003I ). 



We assume that we will be able to split our LRGs into subsamples such that for each 
of them we can def i ne an em pirically calibra t ed te mplate (using a technique similar to that 
of iBudavari et al. I (J2000l ) or iBenitez et al. I (120041 )) which correctly represer its the average 
galaxy colors for that galaxy subsample. Using standard photo-z techniques ( jBenitez Il2000l ) 
we reasonably expect to be able to determine the redshift and spectral type of our galaxies 
in a prelimin ary pass to within . 01(1 + z). This is already being done for the ALHAM- 
BRA survey (IBenitez et al. II2008I : iMoles et al. II2008I ). Thus, for each galaxy we will have a 
preliminary estimate of its redshift to within a^ ~ 0.01(1 + z) and its absolute magnitude 
to with crj\/ ~ 0.15 (within the redshift interval z < 0.9). This ensures that we can pin 
down the required template for each galaxy with large certainty (the error in the absolute 
magnitude corresponds to an intrinsic color variation of only 0.004 mag, much smaller than 
the expected scatter around the sequence at each redshift, 0.03-0.04 mag). 

The LRG template corresponds to z < 0.5 galaxies, and one may wonder if the results 
obtained with this template are representative of higher redshift LRGs. Homeier et al. (2006) 
have measured the V606~-^8i4 colors of a pair of clusters at z = 0.9 with the Advanced Camera 
for Surveys aboard HST. They measure a red sequence that is only 0.09 bluer than the colors 
predicted by our LRG template which illustrates the small amount of color evolution expected 
to z < 1 and shows that the results obtained with our LRG template should be similar to 
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those obtained with real LRG templates at higher redshifts. 

In our simulation we generate the galaxy colors using a combination of the average 
Eisenstein et al. 2003 template, corrected as mentioned above and the first PGA compo- 
nent, multiplied by a coefficient with a Gaussian distribution of rms=1.8. This first PGA 
component scatter represents well the real broadband scatter of galaxies observed by Gool 
et al. (2006) for L > 2.2L^ LRGs. We extrapolate this scatter to L > L^. We also add a 2% 
noise to represent the expected scatter in the zero point determination across the survey 

Then we calculate photometric redshifts using the Bayesian photometric redshift method 
implemented in the BPZ code and a single LRG reference template. We have also tried a 
template library with 11 templates, formed by linear combinations of the LRG template 
and the first PGA component encompassing ±3cr variations and the results are basically the 
same. For simplicity we quote the results obtained with only one LRG template. 

With a single template, there is no point in using a prior, but the Bayesian framework 
still remains useful: it produces the so called "odds" parameter, a highly reliable quality 
indicator for the redshift estimate. In Fig. [TT]we plot the scatter diagram corresponding to 
the quantity {zphot — Zs)/{1 + Zs), where Zg is the true redshift, as a function of the odds 
parameter, together with the rms corresponding to each value of the odds. We can see that 
if we exclude the objects with low values of the odds parameter we get rid of most of the 
redshift outliers. The effectiveness of this techiiique h as been often validated with real data 
( iBenitez II2OOOI : iBenitez et al. 112004 iGoe et al. 1120061). Note that u sing a cut in x^ does not 



work well to eliminate outliers, as it was shown by IBenitez I (J2000l ). 



We thus proceed to eliminate the objects with odds < 0.55 from our catalog. Figure 
shows the scatter diagram for {zphot — Zs)/{1 + Zs) now as a function of the real redshift, Zg. 
Once the odds cut is applied there are no large outliers. 

The resulting redshift and number counts distributions are plotted as solid lines in Figs. [9] 
and [TOl In Fig. [13] we plot the resulting accuracy as a function of redshift. We are safely 
below the 0.003(1 + z) limit for all our redshift range. Finally, in Fig. [THwe plot the number 
density of all the galaxies, and of those with high-quality photo-z as function of redshift. 
These figures show that we have a spatial density oi n > 10"'^ (h/Mpc)^ in the redshift 



range z < 0.9. Since P{k) > 10^ (Mpc/h)^ for LRGs (see eg. Fig. 4 in iTegmark et al. 



(120061 )) and k < 0.2 h/Mpc, we will have nP{k) > 10 for the k range of interest for BAG, 



so that, according to Eq. 1^, shot noise will be negligible. 

Finally, there are two caveats to consider. First, there are no spectroscopic data with 
good enough spectrophotometric calibration for LRGs in the redshift range of interest. We 
can therefore only estimate the intrinsic variation of the galaxies from the data available. 
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We have assumed that it will behave similarly to the variation among LRG types described 
by Eisenstein et al. (2003). Second, the PCA study only covers the 3650 A — 7000 A range, 
and we assume that there is no template variation outside this range. We feel that this is 
justified since most of the redshift information for the galaxies is in practice contained in 
this interval, especially at high redshift. 



3.5. Comparison with a spectroscopic survey 

A typical multi-fiber spectroscopic survey with about 1000 fibers and a resolution R ~ 
2000 in a telescope similar to the one we are assuming here (2-meter class, abo ut 6-deg^ FoV 



etc.) will reach up to a magnitude i < 20 in about 2- hour- long exposures (JBOSS II2008I ). 
assuming the transmission of a good optical spectrograph and low readout noise. This allows 
covering in a year close to 4000 deg^ with 0.1 < 2; < 0.8 for LRGs, or about 2.5 (Gpc/h)^ 
per year. In our PAU approach, with our 300-900 s (depending on the band) exposures, we 
can cover about 2000 deg^ per year with 0.1 < 2; < 0.9 for LRGs, which translates to about 
2 (Gpc/h)^ per year, however with higher galaxy density. This results i n nP(k ) > 10 at the 



relevant scales (see Eq. Q), while for a spectroscopic survey similar to lBOSS I (120081 ) . with 
about 1000 fibers in a 6 deg^ FoV, one can only reach nP{k) ~ 1. Putting volume per year 
and galaxy density together, for an equal-time survey one gets 



(^"/ ")pAU-BAO 




0.6 (22) 



(AP/P)spect 

For the radial modes, one further needs to take into account the slight degradation in infor- 
mation that affects the PAU measurement with its a{z) = 0.003(1 -|- z). 

Furthermore, in the imaging survey one gets many more galaxies than the LRGs. A 
preliminary study for the whole galaxy population obtains a good photometric redshift de- 
termination, a{z) ~ 0.01(1 + z), for a large number of them (over 200 million). These 
galaxies would deliver a constraint on the BAO scale of similar power than the one from 
LRGs (although correlated, since both galaxy distributions trace the same underlying den- 
sity fluctuations), so that the combination of both would improve the sensitivity, and could 
serve as a cross-check on systematic errors. 



3.6. Calibration Requirements 

We present here some general considerations to give an idea of what level of photometric 
and spectroscopic calibration is required to measure the BAO scale with PAU. In the next 
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section we will address the issue of whether these requirements can be met in practice. We 
split this section into photometric and photo-z requirements. 



3.6.1. Photometric Calibration 

The magnitude of a galaxy that we measure in the survey, mo, is the sum of the true 
magnitude m, plus a random statistical error that arises from photon and detector noise, 
emr-i plus a systematic error Cms- The systematic error arises from a variety of effects. For 
example, variations across the survey of the exposure time, mean atmospheric absorption 
and sky background; non-uniformity of galactic dust absorption and inaccuracies in its cor- 
rection; variations in the instrument/detector efficiencies through the duration of the survey. 
All these effects are assumed to have been corrected for through calibrations with standard 
stars and flat fielding corrections. But inevitably, every correction has an error which con- 
tributes to Cms- While the random statistical errors of any two galaxies are uncorrelated, the 
systematic errors in the magnitude over the survey have a correlation function C,ms{z,0,g) 
that is likely to depend on redshift and angular separation 6, as well as galaxy properties g 
(e.g., luminosity, morphology or color). 

The random statistical errors have an effect that is reduced as the number of galaxies is 
increased. Generally the galaxy number shot noise will be larger than the error introduced 
by random errors in the apparent magnitude, as long as these magnitude errors are not very 
large. So if the number of galaxies that is observed is large enough the statistical errors should 
be small compared to the uncertainty in the correlation function due to cosmic variance. The 
systematic errors, however, do not go down with the number of galaxies observed. 

For a flux limited survey, a magnitude calibration covariance across the sky Am{6) 
will result in angular density fluctuations S{6). If we take the number of galaxies brighter 
than magnitude m to be A^(< m) ~ iQ^m (^ typically a In 10 ~ 1), then a magnitude error 
translates into a number density fluctuation error: 

5-a InlOArr, . (23) 

We can decompose the calibration error fleld in the sky into spherical harmonics. We would 
like the resulting spectrum of calibration errors C/" 

C[" = 2n f dcos9A„,{9) Pi {cos 9) (24) 

to produce errors in the angular power spectrum C; which are smaller than the sampling 
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variance errors in C;: 

(jm ^ ^Cl ^ Cl 

' alnlO alnl0y/Uy{l + l/2) 

We will assume that angular clustering will be sampling variance rather than shot-noise 
variance dominated. We also assume Gaussian statistics. The corresponding errors in the 
correlation function Am{6) are: 

2/ -I- 1 

^m{e) = J2^crPi{cose) (26) 

The BAO scale projects at angles between 3.7 and 1.7 degrees for redshifts between 
z = 0.4 and z = 1.0 (smaller redshifts cover a negligible volume). Unfortunately the field of 
view (FoV) of the planed PAU camera plans to cover very similar angular scales. We therefore 
need to be careful about calibration on the FoV. At angular scales 3.7 and 1.7 degrees, the 
requirement in Eq.( l26l) translates into rms correlated calibration errors smaller than 2% to 
3% (ie 0.02 and 0.03 rms magnitude errors) in units of {h/2)/{a In 10) for fsky = 0.2. This is 
for the whole (flux limited) sample (mean z ~ 0.7). These constraints become looser when 
we split the sample into redshift bins because the amplitude of clustering increases as we 
reduce the projected volume. The detailed constraints are shown in Fig. [151 



3.6.2. Selection effects on ^2{j) 

Another angle to look at possible photometric calibration effects is to assume that 
different systematics on the galaxy density fluctuations will act as multiplicative correction 
over the galaxy density at a given position in the sample. We will assume that this type of 
error is uncorrelated to the galaxy clustering so that: 

Us{r) = ar) + Ur) (27) 

where ^(r) is the true correlation and ^e(^) is the correlation due to systematics in selection 
and calibration. To see how this could affect the BAO scale measurement we model the true 
correlation around the BAO scale as a Gaussian peak of width o"o. We further assume a 
generic power-law ^e('") oc r~^ for the error around the BAO scale. The relative shift in the 
BAO scale can be found by Taylor expansion around the peak: 



ie{rBAo) C^i 



Abao = /5 ;r ,\ " ■ (28) 



^{tbao) r 



2 
BAO 
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This requires an amplitude of ^e{rBAo) < 0.002 if we want a shift in the peak Abao < 1% 
and /5 ~ 2. We have used here ctq ~ 15 Mpc/h, rBAo — 100 Mpc/h, and ^{rsAo) — 0.01 
from FigJ31 This corresponds to a 20% error in the correlation at the BAO scale, and about 
\/l^e) = 4.5% error on density fluctuations. 

We have also tested the above calculations directly in simulations, by adding Ce{r) and 
recovering the BAO scale. 

Using Eq. (12^ . this value of C,e{rBAo) < 0.002 corresponds to Am < 0.05 (in units of 
a In 10 with A^(< m) ~ 10"™). Thus, with very different assumptions we reach the same 
conclusion on the requirement on photometric accuracy of around 5% on the BAO scale. 



3.6.3. Photo- z bias 

Systematic errors in the radial direction (photo-z biases) also need to be under careful 
control. At any given redshift, we would like the mean in the photo-z measurements to differ 
from the true redshift by less than 1% (the target in Abao accuracy) in the radial BAO 
distance, ie 0"^ — IMpc/h, which corresponds to: 

A^ = ar H{z)/c ~ 5 X 10"^ {z = 0.8) (29) 

where the numerical value corresponds to ^ = 0.8 and Qm = 0.2. This is about an order 
of magnitude better than the statistical error at the same redshift, i.e. a^ ~ 0.003(1 + z) ~ 
5 X lO'l 

Note that this is a conservative approach because we need the Abao accuracy as mea- 
sured by galaxy density fluctuations and not by the absolute distances to the galaxies. The 
former will probably result into a weaker constraint for A^. 



3.7. Calibration Plan 

As a summary for the above requirements, we need relative calibration to be better than 
about 3% — 5% to avoid systematic effects on density fluctuations to dominate over intrinsic 
fluctuations on the BAO scale. On top of this we would like to have a bias in the photo-z 
scale to be below 1% on radial measurements of the BAO scale. 

In terms of global photometry, it has now been shown that a homogeneous globa l relative 



calibr ation below 2 — 3% accuracy is possible in current and future surveys (e.g., ISterken 



(120071 ) and references therein) . The large fleld of view required by the PAU survey and the 
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drift scanning strategy will both help in the provision of standard calibration techniques, 
such as done in SDSS. We will also need to use a set of calibrated standard spectra of stars 
(or galaxies) to monitor and correct for relative color bias between narrow bands. 

Apart from these "classical" techniques, it seems to be possible to use the observed 
colors of ga laxies with spectr oscopic redshifts as a photometric calibrator. In HSTs Ultra 
Deep Field ICoe et al. I (120061 ) have been able to calibrate the NICMOS zero-points using 
the comparison between predicted colors using the templates of iBenitez et al. I (120041 ) and 
50 spectroscopic redshifts. Similar techniques have been used for the COS MOS field 



(ICapak et al. 1 120071 ) and are being applied to ALHAMBRA (IMoles et al. Il2008f ). Further 



work with the later survey will help refine our calibration redshift requirements. 

Independently of the exact photo-z method finally used for the survey, it will be equiv- 
alent to defining a function z = f{p,C,o), where p are a set of parameters describing the 
function, o are a set of observables like the approximate redshift of the galaxy (determined 
with standard photo-z techniques), its luminosity, size or environment density, and C are 
the observed colors. The LRG population under study is relatively homogeneous and its 
chang es with redshift and rn agnitude can be described with a very compact set of param- 
eters (JEisenstein et al. II2003I ). If we determine the parameters p with enough precision to 
reach a redshift error a^ over the whole range under consideration, then, provided that the 
parametrization is flexible enough to adapt itself to the observed redshift/color relationship 
of galaxies, the systematic zero point error, averaged over all the galaxies, will be equivalent 
to (Jz/\/N^ where az is the rms redshift error and Nc is the number of calibrators. We expect 
to have a calibration set with several thousand galaxies, easy to reach with 10-m telescopes 
as the GTC, and therefore will have negligible redshift bias. 



4. The Survey Instrument 

The approach of the PAU project is to use known and proven technologies to build a 
large field of view (FoV) camera and mount it on a telescope that is optimized for the survey. 

The simulations presented in previous sections use a baseline concept of camera plus 
telescope with a large etendue, AQ ^ 20 m^deg^. This etendue is achieved by means of a 
telescope with a 2-m effective aperture and a camera with 6 deg^ FoV. A total of 42 filters 
is considered, each one having a width of 100 A in wavelength. The full filter system covers 
a range that goes from ~4000 A to ~8000 A, completed by two broad-band filters similar 
to the SDSS u and z bands. 



In this section, we present the main ideas of a possible implementation of such a system. 
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The goal is to show the feasibihty of the telescope/camera system, and not to present a 
complete design. 



4.1. Optics 

Achieving a many-filter, very-large-area survey in a relatively short time and the need 
for a rather large value of the etendue demand a very large FoV. The depth is not a major 
concern since the targets are bright galaxies in all the surveyed redshift range. Therefore, 
the field of view is the main driver of the optical design. 

To give some quantitative estimates, with a telescope of 2.3 m aperture, a pixel size of 
0.4" and state of the art detectors, it is possible to reach S/N ~ 5 for a star of m(AB) = 
23.5 in ~300 sec in all the spectral range bluer than 7500 A. Since the survey should cover, 
as argued before, an area of at least 8000 square degrees, a FoV of 6 sq. deg. is needed to be 
able to perform the survey in 4-5 years. 

It is important to notice that we do not intend to use detailed information on the 
morphology/ shape of the objects, which relaxes the requirements on the image scale. This 
is the reason to choose a rather modest plate scale that translates 15 /im pixels into 26 
arcsec/mm. 

These FoV and plate scale are the basic requirements for a telescope that will be ded- 
icated to the survey until its completion. These optical requirements for such a large FoV 
telescope and the corresponding panoramic CCD camera are demanding but they appear 
feasible. 

The next optical elements are the filters. They are intended to have transmission curves 
with very sharp limit s and minimal waveleng t h overlap, very sini ilar to the filters in the 



ALHAMBRA survey flMoles et al. IboOSl . I2OO8I : iBenitez et al. Il2008h but 100 A-wide 



The location of the filters in the path to the detectors can affect the final efficiency of 
the system. Two options are being considered: attaching the filters directly over the CCDs 
or on plate holders that could be interchanged. The first option is mechanically simpler but 
reduces the survey fiexibility. The second one allows the optimisation of the exposure times 
using different sets of filters depending on the moon-phase or any other external constraint, 
but its practical implementation is more demanding. The final decision will be taken when 
all the practical aspects of the survey, such as observing mode and calibration strategy are 
fixed. 
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4.2. Focal-Plane and Observing Strategy 

The baseline concept for the camera is a large mosaic of CCDs covering the 6 deg^ 
FoV. The scientific goals can be reached with pixels of 0.40". Since most of the current 
astronomical large CCD detectors have pixel scales of 15 microns, we need a camera of 
around 500 Mega- pixel or a number of 2Kx4K CCDs that ranges in between 60 and 80, 
including a few CCDs for focusing and guiding purposes. 

The baseline CCDs under consideration for the PAU camera are the fully depleted, high- 
resistivit y, 250 micron thick de vices developed by the Lawrence Berkeley National Laboratory 



(LBNL) ( jHoUand et al. II2003I ). These CCDs ensure a very high quantum efficiency in the 
red zone of the wavelength region covered by PAU. However, the different possibilities of 
optimisation have some impact in the focal plane instrumentation. There is the possibility 
of having two different types of CCDs covering different regions of the focal plane, in direct 
correlation with the filters, in order to maximise the sensitivity in the whole wavelength 
range. Thin blue- sensitive CCDs correlated with blue filters and thick red-optimized CCDs 
correlated with red filters. Several suppliers of CCDs are available for thin blue-optimized 
CCDs. A final decision for the focal plane intrumentation will be taken considering the 
global optimisation of the survey. 

The camera vessel will need to contain a liquid nitrogen reservoir to maintain the focal 
plane cold for several hours, ensuring stability and not compromising the efficiency of the 
observations. We will also investigate the choice of cryo-coolers or pulse tube coolers, taking 
into account the mechanical requirements that impose a stable precision of a few microns on 
the positioning of the CCDs in the focal plane. 

The best observing strategy for a project of these characteristics is the time-delay-and- 
integrate (TDI) drift scanning mode. Drift scanning is a powerful imaging technique in 
which the telescope is kept stationary and one lets the sky image drift across the CCDs. 
Normally, TDI mode is operated at sidereal rate, and the lines of the CCD must be read 
in perfect synchronisation with the movement of the sky at the focal plane. In this way, 
long continuous strips of the sky are imaged and large fields can be explored automatically, 
which makes this observing strategy particularly well-suited for large surveys. Therefore, it 
is considered as the baseline strategy for PAU. 



4.3. Front-End Electronics 

The readout of a focal plane of the size of the PAU camera is a challenge. One of 
the most attractive options to read out such a large number of CCDs is to use the open 
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system MONSOON (jStarr et al. II2004I ). developed by the instrumentation division of NOAO 



(National Optics Astronomical Observatory, supported by the US NSF). A custom-made 



system or a commercial one, such as the Leach controller (JLeach et al. 1120001 ). are also being 
considered. 

MONSOON is a generic readout system which consists in three kind of boards: the 
Acquisition Board, responsible for the bias voltage generation and the digitization of signals 
coming from the CCDs, the Clock Board, responsible of the clock signal generation needed 
to read out the CCD, and the Master Control Board, responsible of the event building and 
the data transmission to the DAQ computer. The system can be customised to meet the 
specific and demanding requirements of the PAU camera. 



5. Science Capabilities 

5.1. The BAG scale 

We have performed extensive simulations of the PAU survey, assuming 8000 deg^ covered 
up to z = 0.9 using luminous red galaxies with the galaxy density given in Fig. [THand the 
redshift precision from Fig. [131 This results in measurements of H{z) ■ Vg and dA{z)/rs, Vg 
being here the sound horizon at recombination, which are plotted in Fig. [161 The relative 
precision achieved improves monotonically with increasing redshift, flattening out at about 
5% for H{z) and 2% for dA{z). 

We have split the redshift interval from z = 0.1 to z = 0.9 into 16 equal bins. Results do 
not change if we change the binning. Combining these measurements into a cosmology fit, 
taking into account correlations, leads to determinations of th e parameters r2m,Wn,Wn- In 



the following , we will be using the standard parametrization (jChevallier &: Polarski 112001 



hinder II2003I ) of the time evolution of the Dark Energy equation of state, w{z) = Wo+Wa-{l — 
a), where wq denotes the equation of state now, and Wa is (minus) its current derivative with 
respect to the scale factor a. The value of the reduced Hubble constant, h drops out of the 
measured quantities. The left panel of Fig. [T71 shows the 68% confidence level (CL) contour 
in the flm-w plane that can be achieved using only PAU LRG data. The corresponding 
one-sigma errors are a{flm,w) = (0.016,0.115). A flat universe and constant equation of 
state has been assumed. In the right panel of Fig. [T71 68% CL contours are shown in the 
WQ-Wa plane, assuming a fiat universe. The outermost contour approximates the expected 
world combined precision from SNe and WMAP when PAU will start taking data, while 
the inner contour adds the PAU LRG data to the previous data set. The reduction in area 
corresponds to an improvement by more than a factor three in the DETF figure of merit. 
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The one-sigma errors are a{wQ,Wa) = (0.14,0.67). 

We have also simulated a straw-man spectroscopic survey with equal area and depth 
but with az = 0.0005(1 + z). The greatly improved redshift precision results in only a mod- 
est 20% decrease in the cosmological parameter uncertainties. Actually, the larger galaxy 
density that a photometric survey affords overcompensates for this deficit, as can be seen in 
Fig. [IHl where our simulated reach for several proposed BAO surveys is compared. Details 
of the survey characteristics are taken from public sources and are summarized in Table [2J 
For Wiggle-Z, we adopt the total number of galaxies and redshift distribution as given in 



Glazebrook et al. I (120071 ). Surprisingly, our Wiggle-Z constraints on dark energy parame- 
ters are approximately a factor two worse than what one would predict from the errors on 
the d^ and H measurements quoted by the authors. We use the information in section 3 



of the SDSS3 project description (http://www.sdss3.org/collaboration/description.pdf) to 



simulate BOSS. The information for HETDEX is taken from iHill et al. I (120041 ). The WF 



MOS surveys details are taken from the Feasibility Study Report, which can be found in 



|http://www.gemini. edu/files/docman/science/aspen/ The Pan-STARRS 1 Survey (PSl) 



and Dark Energy Survey (DES) information are taken from their web pages (http://pan-starrs.ifa.hawaii.edi 
and [http: //www. darkenergysurvey.org/ 1 ) . 



5.2. Other Probes of Dark Energy 

5.2.1. The galaxy clustering in redshift space 

The redshift accuracy will be sufficiently good to identify individual structures (walls 
and voids) along the line of sight (see Fig. 2). Accurate measurements of the redshift-space 
power spectrum will be possible in the linear and mildly non-linear regime, and a detailed 
comparison with theoretical predictions will be done in conjunction with the measurement 
of Baryon Acoustic Oscillations. Detailed measurements of redshift space distortions offer 
an independent test of the growth history of the peculiar velocity field (i.e. Newtonian 
gravitational potential). This encloses cosmological information on dark energy (and/or 
modified gravity), complementary to that in BAO, which measures the background history. 
A descomposition of the 2-point correlation function in radial and transverse directions also 
allows for a measurement of bias h and the amplitude of matter clustering as, which can 
be used to study the growth history of density fiuctuations to 2; = 1. Measurement of the 
amplitude of the galaxy power spectrum, P(A;), as a function of redshift can also be used 
to determine the growth rate of structure through the cross-correlation of the galaxy data 
with future CMB lensing data or by using higher-order correlations to determine the bias 
parameter as a function of scale and redshift. Higher-order correlations, such as the 3-point 
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function or bispectrum, can also be used to measure the BAO feature. 



5.2.2. Weak Lensing 

Weak lensing is sensitive to both the distance and the growth factor as a function of 
redshift. The lensing effect can be measured using either the shear or the magnification. The 
PAU camera will not be optimized to measure galaxy ellipticities, so weak lensing shear may 
not be as good as those from other surveys. However, the accurate photometric redshifts 
obtained in PAU may be combined with ellipticity measurements obtained in other surveys 
for the same galaxies. This additional information will help separating shear lensing from 
intrinsic galaxy aligment. 

Gravitational lensing modulates the observed spatial distribution of galaxies. Dim galax- 
ies that otherwise would not have been detected are brought into the sample by the lensing 
magnification. This increases the observed number density of galaxies. On the other hand, 
magnification also increases the apparent area, which leads to a drop in the observed num- 
ber density of galaxies. The net lensing effect, known as magnification bias, is controlled 
by the slope of the number counts. The PAU survey will be able to measure this effect by 
cross-correlating galaxy samples defined by separated redshift slices. 



5.2.3. Galaxy Clusters 

Galaxy clusters are the largest collapsed structures in the Universe, containing up to 
hundreds or thousands of individual galaxies. The redshift distribution and the evolution 
of clustering of massive clusters of galaxies can provide a direct measurement of the cosmic 
volume as a function of redshift as well as the growth rate of density perturbations. This is 
complementary to the measurement of the BAO scale, which is purely geometrical in nature. 
Comparison of theory to observations requires a calibration of the cluster masses (or at least 
the mass threshold of each cluster sample). Clusters of galaxies can be identified optically 
by searching for concentrations of galaxies with the same color. The PAU survey by itself 
will provide a new window for accurate optical cluster detection and selection, based on the 
combination of photometric colors and good photo-z precision over all galaxies around each 
cluster, which will help improving on cluster completeness and contamination. PAU will 
also provide the opportunity to self-calibrate the mass threshold of a given cluster sample in 
different ways, such as stacking weak lensing magnification measurements over the cluster 
positions or using the (biased) amplitude of clustering in the same cluster sample. The 
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photo-z accuracy for clusters will be improved in comparison to the galaxy photo-z by the 
square root of the number of galaxies in the cluster. This will result in a typical photo-z 
accuracy which is a few times smaller than that for galaxies. At the same time, one could 
use the velocity dispersion of the galaxies in each cluster to provide an estimate of the cluster 
mass. This will not be accurate for individual clusters, but should be accurate enough to 
have an estimate of the mass threshold of a given cluster sample. 

A cluster survey carried out over the PAU area also constrains cosmology through the 
spatial clustering of the galaxy clusters. As mentioned above, this can be done with even 
higher photo-z accuracy than in the PAU galaxy survey. The clustering of galaxy clusters 
reflect the underlying clustering in the dark matter; these correlations contain a wealth 
of cosmological information, much like the information contained in the LRG correlation 
function, including the BAO position. Even if the number density for clusters is lower than 
that of LRGs, this is partially compensated by the higher (biased) clustering amplitude. 
We plan to use the PAU cluster redshift distribution and the cluster power spectrum and 
clustering as cosmological probes to study the density and nature of the dark energy. 

PAU can also be used in combination with other surveys to provide accurate photo-z in 
a sample of clusters detected by the Sunyaev-Zel'dovich (SZ) or X-ray signatures of hot gas 
in clusters, as well as in weak lensing cluster selection. 



5.3. Other science 

The large number of narrow filters in the PAU survey will yield many colors for all 
the detected galaxies, allowing the measurement of numerous interesting parameters for the 
study of galaxy evolution: stellar mass, stellar age distribution, metallicity, dust absorption, 
and interstellar gas emission. This will make possible a detailed study of the rates of star 
formation, galaxy mergers and chemical evolution that can account for the evolution of the 
stellar contents of galaxies of different types, as a function of their environment. 

The PAU survey will substantially increase the sensitivity of astronomical observations 
to the presence of intergalactic dust, and possibly detect and characterize it for the first 
time. Intergalactic dust extinction can be searched for by correlating the foreground density 
of galaxies with the background surface brightness of sources, as well as the extinction 
measured from our multiple colors. Our accurate photometric redshifts will allow for a 
good estimate of the mass column density (which should presumably correlate with dust 
extinction) in the line of sight to every background galaxy. A detailed extinction curve for 
intergalactic dust may be measurable with our multiple narrow-band colors. The presence 
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of gray dust (causing extinction with no reddening) could also be explored. 

The narrow filters of PAU will result in an improved separation of quasars and stars 
compared to other surveys, and also a more accurate estimate of quasar redshifts from the 
photometry. Our narrow-band photometry may also be useful to study the mean transmis- 
sion of the Lyman alpha forest on large scales, its fluctuations and its evolution with redshift. 
Quasar lens candidates can be searched for by checking for the presence of multiple images 
and foUowed-up at another observatory with higher angular resolution. The PAU survey 
will also provide a map of galaxies in the vicinity of the lines of sight of all our quasars, 
with rather accurate photometric redshifts for all the galaxies. One can follow up with spec- 
troscopic observations of a selected sample of quasars at another observatory, and use the 
PAU galaxies to correlate absorption systems seen in the quasar spectra with the galaxies, 
to study the distribution of gas around galaxies. 

The most luminous star-forming galaxies at redshifts z > 2.5 will be detectable with 
the PAU survey narrow-band photometry by means of their Lyman continuum break, the 
Lyman alpha forest absorption, and possibly a Lyman alpha emission line. This will allow 
the study of this galaxy population and its clustering properties over an unprecedentedly 
large volume. 

By selecting the PAU filters appropriately, several parameters of the stars observed in our 
survey should be measurable, such as effective temperature, surface gravity, iron abundance, 
and a/Fe. An accurate determination of the density profile and metallicity distribution of 
halo stars in the Milky Way may follow: the PAU survey could yield the largest number of 
metallicity measurements of halo stars, characterizing the stellar populations of the various 
streams believed to have originated our stellar halo. Giant stars may be also detected by 
PAU at distances up to ~ 1 Mpc. This could provide our first substantial sample of stars 
far from any galaxy in the Local Group, and extend the measurements of the halo profile of 
the Milky Way and its metallicity distributions to much larger radii. For this purpose, it is 
necessary to have good ways of distinguishing nearby K and M dwarfs from distant K and 
M giants with the PAU narrow-band photometry. 

PAU should also be great for serendipitous discoveries. We will have spectral information 
for every one of the 10^ pixels in the survey, which will allow for search of diffuse, low signal- 
to-noise, components or rare new objects. For example, if there exist any objects in the 
universe that produce bursts and emit most of their power in an emission line, the PAU 
survey will be ideal for discovering them. 
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6. Summary and Conclusions 

In 1998, the discovery of tlie accelerated expansion of the universe changed completely 
our understanding of the universe and its components. Ten years on, the quest to understand 
what causes the acceleration continues. Along the way, the Baryon Acoustic Oscillation 
(BAO) technique has been identified as a systematically robust, yet statistically powerful, 
probe of dark energy properties. In particular, measuring the BAO feature along the line 
of sight as a funcion of redshift provides a direct determination of H{z), which, in turns, 
depends on the amount and characteristics of dark energy. However, such a measurement 
requires a very precise determination of galaxy redshifts. 

We have presented here a novel approach to photometric redshift determination which 
allows the measurement of the BAO feature along the line of sight in an efficient way, by 
using a set of about 40 narrow-band (FWHM^ 100 A) filters. The approach complements 
(for BAO and for other science) spectroscopic surveys, which typically measure much more 
precisely the spectra of a much reduced sample of galaxies. 

Because of the intrinsic width of the peak in the galaxy-galaxy correlation function of 
about 15 Mpc/h, there is a fundamental limitation as to how much one can improve the BAO 
measurement by reducing the photo-z errors. A redshift precision of order a{z) ^ 0.003(l+z), 
corresponding to 15 Mpc/h along the line of sight ai z = 0.5, is about optimal for this 
measurement. Redshift space distortions, biasing and non-linear effects produce distortions 
that can be comparable to the effect of this photo-z precision. 

Simulations show that both the target galaxy density (n ~ 10~^(h/Mpc)^) and precision 
in redshift (cr^/(l + 2;) ~ 0.003) can be achieved with the proposed system. These simulations 
indicate that PAU by itself can determine the equation of state of the dark energy assumed 
constant {w) to about 5%, while when the PAU data is combined with expected supernova 
and CMB data samples, a sizable increase in the Dark Energy Task Force figure-of-merit 
(inverse of the area of the error ellipse in the WQ-Wa plane) by about a factor three is 
achieved, making the PAU very competitive when compared to other planned ground-based 
BAO surveys, photometric or spectroscopic. 

The survey will produce a unique data set with low-resolution spectroscopy in the optical 
wavelengths for all objects in the northern sky up to itlj^b = 23-23.5 arcsec"^ (five sigma). 
A survey like PAU, producing such a catalogue, will have enormous legacy value and will be 
extremely useful for many areas of astrophysics, with contibutions that are different from, 
and complementary to, those a spectroscopic survey can deliver. 
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Table 1. N-body simulations used in this paper. Minimum halo mass and number of halos 

correspond to z = 0.5. 



name 



iV„ 



halo mass 



N, 



halo 



acronym Mpc/h number lO^^Msun/h Total number 



MICE3072 
MICE1536 



3072 
1536 



10243 
10243 



>375 

>47 



1.1 X 106 
2.1 X 10^ 



Table 2. Details of the surveys considered in fig. [181 



survey 


z range 


Number galaxies 


Tracer 


Area 
deg2 


Volume 
(Gpc/h)3 


Radial 
information 


Time scale 


reference 


WiggleZ 


0.3< 2< 1.2 


2.8 X 10^ 


ELG 


1000 


2.04 


yes 


2007-2009 


Glazebrook et al. (2007) 


BOSS-LRG 


0.2 < z < 0.8 


1.5 X ID'S 


LRG 


10000 


8.06 


yes 


2009-2014 


see text 


HETDEX 


1.8 < z < 3.7 
0.5 < z< 1.3 


1.0 X 10" 
2.0 X 10" 


LAE 
ELG 


200 
2000 


1.91 

4.47 


yes 
yes 


? 
? 


HiUetal. (2004) 


WFMOS-ELG 


see text 


WFMOS-LBG 


2.3 < z<3.3 


6.0 X 10'""' 


LBG 


300 


1.53 


yes 


? 


see text 


PSl 


0.3< z< 1.5 


5.0 X 10* 


ALL 


20000 


65.3 


no 


? 




DES 


0.3< z< 1.5 


1.5 X 10* 


ALL 


5000 


16.3 


no 


2011-2015 




PAU-LRG 


0.1 < z<0.9 


1.3 X lO'^ 


LRG 


8000 


8.6 


yes 


2011-2015 


this paper 
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Fig. 1. — Relation between the change in the dark energy equation of state parameter w, 
shown in the x-axis, and its effect in the apparent measured BAO scale, Abao^ shown in the 
y-axis. Changes are shown in percent, relative to the w=-l case. Top panel corresponds to 

AH{z)/H{z). The bottom panel shows the angular diameter 
AdA{z)/dA{z). The different lines correspond to z = 0.3 (continuous), 
z = 0.6 (short dashed) and z = 1 (long dashed). In all cases Qm = 0.25 and fiat universe are 
assumed. All other cosmological parameters are kept fixed. 



the radial distance: A^^q 
distance: A^^q 
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Fig. 2. — Sytematic effects in the lightcone: panels show a 1 Mpc/h thick section of the 
hghtcone distribution in MICE3072 in comoving coordinates. The two bottom panels corre- 
sponds to the actual dark matter distribution in the simulation in real (bottom) and redshift 
space (upper panels). The top two panels also include a (Gaussian distributed) photo-z error 
distortion of cr^ = 0.003(1 + z), as expected from PAU galaxies, and an order-of-magnitude 
worse case, a^ = 0.03(1 + z). The BAO scale is shown by a section of circle with radius 100 
Mpc/h around the observer. 
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Fig. 3. — Smearing of the BAO signature due to photometric redshift errors. The cir- 
cles denote the two-point correlation function from over a million halos with mass M > 
3.7 X 10^^ h~^ Mq (assumed to host LRGs) measured in a MICE simulation of 27 h~^ Gpc^ 
volume. The dashed line is the linear correlation function scaled with the linear halo bias 
(b = 3), while the black solid line corresponds to the nonlinear prediction given by RPT 
( Crocce fc Scoccimarro II2008I ). Their difference shows the degradation coming solely from 
nonlinear clustering. In addition, the triangle (red), square (blue) and cross (green) symbols 
show the measured correlation function after a Gaussian error degradation in the line-of-sight 
position of the halos is introduced (crz/(l + z) = 0.003,0.007 and 0.03, respectively). The 
corresponding solid lines are the analytical predictions derived from Eq. ([8]). The right panel 
shows a zoom over the peak region scaled as r'^C,{r). Clearly, the signal-to-noise in the BAO 
feature reduces with photo-z error and starts to totally disappear above the PAU threshold 
of 0.003, which roughly corresponds to the intrinsic width of the BAO feature due to Silk 
damping. 
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Fig. 4. — Ratio of the power spectrum measured at z = 0.5 in MICE1536 to a smoothed 
version of the same spectrum for dark matter (left panel) and halos of mass M > 4.7 x 
10^"^ h^^ Mq (right panel). These large halos are expected to host LRGs. The reference 
spectrum was obtained by rebinning the original one in bins of Ak = 0.055 /iMpc~^ in 
order to wash out the BAO signature but keeping the broad band shape of the nonlinear 
spectrum. Error bars were obtained using the approximation in Eq. (Q. The red solid line 
corresponds to the parametric fit given by Eq. (ITOl) with A/vbao = 0.016 for dark matter 
and A/tbao = 0.017 for halos {tbao = 108.6 Mpc/h for our reference cosmology). This 
figure illustrates that for both dark matter and halos, one can approximately model BAO 
in the P{k) with E g. (ITOl). This conclu sion also applies in redshift space and for different 
galaxy populations (JAngulo et al. II2008I ). 
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Fig. 5. — An example of a filter system similar to the one which will be used by the PAU 
survey. We have included the redshifted spectrum of an early type galaxy at z=0.2 from the 
Bruzual and Chariot library to illustrate how the sharp 4000 A break (which here falls at 
4800 A) is basically bracketed by only two filters. Note that the filters are spaced 93 A but 
have FWHM widths of 118 A due to the wavelength extent of their wings. The blue squares 
represent the flux which would be observed through the filters. Note that many spectral 
features apart from the 4000 A break are resolved by such a filter system. 
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Fig. 6. — The sky background assumed for our simulations. We have assumed that we 
are able to adapt the choice of filters in our observations to the moon cycle, observing in 
the u-band the darkest night, and then moving towards redder filters as the sky brightness 
grows. The red, continuous line corresponds to our expected "effective" sky spectrum, with 
the squares showing the equivalent broad band AB magnitudes in the UBVRI filters . The 
spectrum is normalized to have the same broad band brightness as the iPatat I (J2004j ) mea- 
surements of the dark sky at Paranal for the U and B bands, and the same as the middle of 
the cycle "gray" nights from IWalker I (119871 ) in the rest of the filters. 
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Fig. 7. — Expected limiting magnitudes for PAU-BAO. The squares represent the 5a mag- 
nitude hmits within a 1 sq. arcsec aperture, the circles within a 2 sq. arcsec aperture and 
the continuous line is the 5a magnitude limit which would be reached if we divided the total 
exposure time of 19440s equally among all the filters. 
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Fig. 8. — The expected 5-0" limiting magnitudes for point sources (squares) and the observed 
spectra of a L^ red galaxy at different redshifts (without taking into account spectral evolu- 
tion, but taking into account aperture corrections). Note that we are able to catch the rest 
frame 4000 A break with enough filters on both sides up to z=0.9. 
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Fig. 9. — Differential number counts distribution of L > L^ red galaxies. 
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Fig. 10. — Redshift distribution oi L > L^ red galaxies. 
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Fig. 11. — Photometric redshift error as a function of the Bayesian odds. Note that a cut at 
odds = 0.55 ehminates most of the objects with high redshift errors. For the sake of clarity, 
only one in every five points is plotted. The solid line corresponds to the rms of ^z/{l + z) 
for each value of the odds. 
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Fig. 12. — Scatter plot comparing the normalized difference between the photometric red- 
shifts and the "true" input redshifts zs- The red points are eliminated by the odds < 0.55 
cut. For the sake of clarity, only one in every five points is plotted. 
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Fig. 13. — Photometric redshift error as a function of redsliift, for all L > L^^,I < 23 red 
galaxies, and for tlie subset witli liigii quality photo-z. 
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Fig. 14. — Spatial density as a function of redsliift, for all L > L^,,I < 23 red galaxies, and 
for the subset with high quality photo-z. 
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Fig. 15. — Required systematic calibration error (rms percentage) for two broad redshift 
slices (thick blue continuous lines): z = 0.4 — 0.5 (top) and z = 0.9 — 1.0 (bottom) and for 
a flux limited sample (red dashed line) including all galaxies to the depth of PAU (mean 
z = 0.7). In these units, at BAO scales (which is a function of z and is marked by the arrow) 
the correlation in calibration error has to be smaller than about 6% for z ~ 0.45 and 5% 
for z ~ 0.95. For other science, the stronger requirements are driven by the flux limited 
sample, i.e. < 2% and < 8% in correlated errors on scales smaller than 4 and 0.1 degrees 
respectively, as given by the dashed line. 
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Fig. 16. — Left: The expected measurement of radial BAO scale from the PAU survey 
(LRGs only). Right: Same for the measurement of the transverse (angular) BAO scale. 
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Fig. 17. — Left: 68% confidence- level contours in the flm-w plane, using only PAU data, 
assuming a fiat universe and a constant equation of state w. Right: 68% confidence-level 
contours in the WQ-Wa plane for the world combined data from SNe and WMAP in about 
2010, and after adding PAU data to that data set. The area of the contour decreases by 
about a factor three. A fiat universe has been assumed. 
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Fig. 18. — Precision on w (assumed constant) for different proposed ground-based BAO 
surveys. All other cosmological parameters are kept fixed, therefore the overall scale is 
unrealistic, but the relative reach of the different proposals should be realistic. Details 
about the inputs for the calculations can be found in Table O 



